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Abstract 

In the previous paper hep-th/0604112 we calculated the first of the five planar two- 
loop diagrams for the Lcc vertex of the general non-Abelian Yang-Mills theory, the vertex 
which allows us in principle to obtain all other vertices via the Slavnov- Taylor identity. The 
integrand of this first diagram has a simple Lorentz structure. In this letter we present the 
result for the second diagram, whose integrand has a complicated Lorentz structure. The 
calculation is performed in the D-dimensional Euclidean position space. We initially perform 
one of the two integrations in the position space and then reduce the Lorentz structure to D- 
dimensional scalar single integrals. Some of the latter are then calculated by the uniqueness 
method, others by the Gegenbauer polynomial technique. The result is independent of the 
ultraviolet and the infrared scale. It is expressed in terms of the squares of spacetime intervals 
between points of the effective fields in the position space - it includes simple powers of these 
intervals, as well as logarithms and polylogarithms thereof, with some of the latter appearing 
within the Davydychev integral J(l, 1, 1). Concerning the rest of diagrams, we present the 
result for the contributions correponding to third, fourth and fifth diagrams without giving 
the details of calculation. The full result for the Lcc correlator of the effective action at the 
planar two-loop level is written explicitly for maximally supersymmetric Yang-Mills theory. 
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1 Introduction 



It has been shown in Refs. [1, 2, 3] that the effective action of dressed mean fields for J\f = 4 
super- Yang-Mills theory does not depend on any scale, ultraviolet or infrared. These results 
were derived from the results of Refs. [4, 5, 6, 7, 8]. Scale independence suggests that kernels 
of these dressed mean fields can be analyzed by the methods of conformal field theory. We 
started to investigate the simplest scalar vertex Lcc in the Landau gauge in Ref. [9], where c 
is a real ghost field and L is the auxiliary field which couples at the tree level to the BRST 
[10, 11] transformation of the c field. This vertex is simple in the Landau gauge where it 
is totally finite for M = 4 super- Yang-Mills theory. By solving Slavnov- Taylor (ST) identity 
[12, 13, 14, 15, 16, 17], all other vertices in that theory can be derived from this vertex [4]. The 
ST identity is a consequence of the BRST symmetry of the classical action [10, 11]. Recently, 
by using unitarity methods it has been demonstrated up to the four-loop level that the only off- 
shell conformal integrals in the momentum space contribute in the maximally-helicity-violating 
(MHV) four-particle amplitudes [18], and iterative structure has been conjectured for all the 
MHV amplitudes [19]. 

As has been found in Ref. [20], all the contributions to the off-shell four-point correlator 
of gluons in that theory up to three-loop level (at least) are equivalent to ladder integrals (UD 
integrals) calculated in the momentum space by Usyukina and Davydychev in Refs. [21, 22]. 
They are expressed in terms of certain functions (UD functions). It is known from Refs. [21, 22, 
23, 20] that UD functions are conformally invariant in the momentum space. Fourier transform 
of UD functions in the momentum space can be expressed in terms of the UD functions in the 
position space which possess the same property of the conformal invariance but this time with 
the arguments of the position space. On the other hand, the scale independence of the effective 
action of the dressed mean fields is a consequence of the vanishing of the beta function in the 
maximally supersymmetric Yang-Mills theory and the vanishing in turn is a consequence of the 
conformal symmetry of the theory. All this suggests that correlators of the dressed mean fields 
can be analyzed by the technique of the conformal field theory. 

To solve the ST identity and to check the conformal invariance of the effective action of the 
dressed mean fields explicitly, we should work in the position space. The two-loop contribution 
contains five diagrams. We calculated the first contribution in the previous paper [9]. Now we 
calculate the second contribution that corresponds to the diagram (b) of Ref. [9] , give the result 
for the rest of the diagrams, and present the full result for the planar two-loop contribution to 
the Lcc correlator of the effective action in the Landau gauge for Af = 4 supersymmetric Yang- 
Mills theory. The notation used here is the same as in Ref. [9]. The Lcc vertex is superficially 
convergent in the Landau gauge. This fact can be checked by index counting and by noting 
that two derivatives from the ghost propagators can always be integrated out of the diagram 
due to the transversality of the gauge propagator. It means that the field c does not have 
renormalization in the Landau gauge. Formally, this result holds to all orders of perturbation 
theory due to the so called antighost equation [24] . In the nonsupersymmetric theories this vertex 
is not finite and a calculation of the anomalous dimension of operator cc has been performed in 
[25, 26]. Finiteness of this Lcc vertex at planar two loop level in Landau gauge in maximally 
supersymmetric Yang-Mills theory is equivalent to the cancellation of poles between fourth and 
fifth diagrams of Ref. [9] that is, in turn is equivalent to the vanishing of the (3 function of the 



1 



gauge coupling in that theory. 

Knowing the structure of the Lcc vertex one can obtain the structure of other irreducible 
vertices by solving the ST identity. The algorithm can be applied also to other theories different 
from the theory in consideration. Among possible applications are M = 8 supergravity, Chern- 
Simons theory, string field theory, massless gauge theory near fixed points in the coupling space, 
and topological field theories in higher dimensions. 

The paper has the following structure. In Sec. 2 we write the integral expression that 
corresponds to the diagram (b), and analyze it by dividing it in three parts. In Sec. 3, the result 
for the first part is reduced to single D-dimensional integrals (D = 4 — 2e). The analogous 
expressions for the second and the third part are written in Appendix A. The sum of all three 
parts contains single D-dimensional integrals J(a±, ct2, as) with the sum of indices a± + «2 + «3 
in the denominators equal to D — 1, D, D + 1, D + 2. In Sec. 4 we present the part with the 
integrals whose sum of indices is equal to D — 1. In Appendix B we present the part with the 
integrals whose sum of the indices is D, D + 1, and D + 2; these integrals are calculated there 
explicitly by using the uniqueness method and its variants [27, 28]. In Sec. 4 we reduce the 
aforementioned {D — l)-type partial sum to explicitly known terms and to terms containing 
the integrals J(l, 1, 1), J(e, 2 — 3e, 1) and derivatives thereof - by applying formulas obtained 
in Appendix C where the method of integration by parts procedure (IBP) [29, 30] was used. 
In Sec. 5 we obtain the entire result for the diagram (b) in D dimensions, in known terms 
and terms with the integrals J(l, 1, 1), J(e, 2 — 3e, 1) and derivatives thereof. The terms with 
J(e, 2 — 3e, 1) are calculated explicitly in Sec. 5, while derivatives of J(e, 2 — 3e, 1) are calculated 
in Appendix D - in both cases those formulas of Appendix C were employed which were obtained 
by the Gegenbauer polynomial technique (GPT) introduced originally in Refs. [31, 32, 33, 34] 
and further developed in Ref. [35]. In Appendix E, all the contributions for the diagram (b) in 
D dimensions are collected, and the e — > limit is performed. In Sec. 6 we present this result 
in a more explicit and shorter form. In Sec. 7 we write the result of calculation of the rest 
of diagram and present the total result for the planar two-loop Lcc correlator of the effective 
action in the Landau gauge. In Conclusions, we comment on the result obtained and on further 
possible developments of the analysis performed in this paper. 

2 Diagram (6) 

The two-loop correction to the Lcc vertex can be represented as a sum of five diagrams depicted 
in Fig. 1. We have introduced for brevity the notation 

[yz] = (y - z) 2 , [yl] = (y-xi) 2 ,.... 

and so on. The structure of Lorentz indices in diagram (b) is more complicated than the structure 
of the indices in the diagram (a) considered in Ref. [9]. The algebra can be inferred immediately 
from the diagram (6). 

(gg)g ( 9aX , 9n _ A2y) a (2y) x \ (yl) x (zy)^ f _ (3z)^3z) u \ (31), 

[2z]*-< W" e N 2 - e / [yi] 2 - e N] 2 - £ W" e+ 1 ; [M 2 ~ e y [3i] 2 ~ e 

We work in the position space, in D = 4 — 2e dimensions, using the technique of dimensional 
regularization [36], thus maintaining the intermediate results finite. We have to calculate the 
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Figure 1: The two- loop diagrams for the Lcc vertex. The wavy lines represent the gluons, the straight lines the 
ghosts. Black disk stands for the total one-loop correction to the gluon propagator. 



2D-dimensional integral 1 

(2z) a / g aX _ (2y) tT (2y) x \ (yl) x 

\2z\ 2 - W" £ K } [2y] 2 ^ J [ y l]2- 



I 



(zy)n ( 9i*v , on ^(3z)^(3z) u \ (31) 



+ 2(1 



[zy] 2 -' VN 1 - 6 v ' [3z} 2 -< J [31] 2 - £ 



By using the representation (yl)\ = (y2) x + (21) x , after simple algebra we can represent this 
integral as a sum of three integrals 

(2z) a (2y) a (zy)^ ( g^ u (3z) fl (3z) u \ (31) v 



-(1-0 \ [ D n (2z) a (2y) a (zy)^ / g^ _ i 

1 >J V [ylf-^zf-^yY-^zyf-AlM 1 -" e) '[3z} 2 ^ J [31] 2 - 

( n„ n? [(2^) g (2j/)q][(2yA)(2i)A] {zyh. ( _g EL _ , 9n ^(3^(3z)^ (31), 

+( >J V [yl] 2 ^[2z] 2 ^[2y]^ [zy] 2 ^ \ [3z]^ + ' U 



[3z] 2 ^ J [31] 2 ^ 
(2z) a (21) a (zy)^ ( 9lxv J3z)^3z) u \ (31), 



.[ Dv Dz (2^(21). j^/W^_ + 2(1 



[3z] 2 " e J [31] 2 ^" 



As in the previous paper [9], we modify the e's appearing in the propagators (e i— > Kj-e) in order to 
make the uniqueness method applicable and thus to make the calculations easier. Such changes 
do not affect our result in the e — > limit since the integral is finite. We do the following 
changes: in the powers of the denominators of the four ghost propagators in Eq. (1), we change 
e i— > 2e in the first one, and e i— > in the other three; in the powers of the denominators in the 
two gauge propagators, we change e i— ► 2e in the first and e i— > in the second propagator; in the 
nominator of the second gauge propagator, we change e i— ► e, where e is as in the previous paper; 
in the nominator of the first gauge propagator, we change e i— > I2 where €2 is different because 
it is related to the requirement of keeping the transversality in the modified propagators [9] 

e --(i^)' e2 -(iT^r (2) 

The above integral for the diagram (b) now acquires the following form: 

(2z) a ( g uX , o/i _ . (2y) CT (2y) A \ (yl) A 



We use the D- dimensional measure Da; = 7T-^-d D a; [9]. 
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v (zy)» ( 9^ „ n > (32) M (3z)„ ^ (31)„ 
^^l[3l] 4 



VN v ; [3z] 2 7 [3i] 2 

= -(3-2e 2 )/ 6)1 + (2-2e 2 )/ 6)2 + 7 M , (3) 

where the three integrals I^j are 

f n n (2^)<x(2y) CT (*y)„ / 

v M - y uz [ yl ]2[ 2z ]2-2 e [ 2y ]i-2 e [zy] 2 ^ + ^ - e J 

r y n , p . [(2^(2y) g ][(2y A )(21) A ](^) M /g^ 

r /* n n (2z) CT (21) CT (zy)^ / ^ 

In these 2D-dimensional integrals, we first perform the Lorentz algebra only for certain subprod- 
ucts in the numerators, then perform one of the two D-dimensional integrations by using the 
uniqueness method or its variants, 2 and then finish the Lorentz algebra. We found this method 
simpler than making the entire Lorentz algebra initially and then integrating the contributions. 
After integrating, the number of terms in the algebra is significantly reduced. 



(3z) M (3z)„\ 


(31), 


(4) 


[3^] 2 J 


[31] 2 ' 


(3z)^(3z 
' [3^] 2 


)u\ (31), 
/ [31] 2 


(5) 


(3*)„(3z)„\ 


(31), 


(6) 


[M 2 J 


[31] 2 • 



3 Reduction to single scalar integrals 



Following the afore-described procedure, we can scalarize the expression for each of the three 2D- 
dimensional integrals in Eq. (3), reducing them to linear combinations of integrals J(ai, a 2 , 03) 
where 

J(ai ' a2 ' a3)= /^[xl]-[x2]-[-3]- • 
For the first integral we obtain 

,4(l,l-2e,l) 17 6(l-e 2 )(l-2e) 2(1 - e 2 )(l - 2e) ' 



hi 



4e(l - 2ef 



[31] 2 [12]- 



+ 
+ 



2(l- e 2 )(l-2 e) 
+ [31] 2 [12]i- J[€ ' 2 d£,ij + 
2(l-e 2 )(l-2 £ ) 



[31] 2 [12]!- e 
6(l-e 2 )(l-2e) 



[31] [12]- 



J(2 + e,l -3e,0) 
J(2 + e,2-3e,l) 



[31][12]i-« 




2(l-e 2 )(l- 


-26) 


[31][12]i- 


■e 


(1-26 2 )(1- 


-26) 


[31] 2 [12]i 


-e 


2(l-e 2 )(l- 


-26) 



+ 



e — e 



[31] 2 [12]!- e 

€ — € 



I=7 J(2,l-3e,l + e) + 



[31] 2 [23] 



J(l + e,2- 


3e,i; 


J(2 + e,l- 


3e,l) 


J(l + e,l- 


3e,l) 


J(l + e,2- 


3e,0; 








3e, e] 



[31] [23] 

4-8e + 3e 2 , N 4 - 8e + 3e 2 , 

-J(2,2-3e,l + e)+ roi12roo1 _ J(l, 2 - 3e, 1 + e 



[31] [23]" 



[31] 2 [23]- 



By the "variants of the uniqueness method" we mean the cases when the sum of the powers in the denominator 
of the integrand is larger than D but can be made equal D once the integrand is represented as a derivative, or 
derivatives, with respect to a coordinate that is not integrated over. 
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/ e (l-2e) 4 - 106 + 66 2 [12]e(2-3e) \ , 
\J31H23] 1 -* [31] 2 [23]-^ [Sl] 2 ^] 1 - 6 J ^ ' ' ; 



+ 
+ 



(l-2e)(l-e) (l-e) 2 + [12](l-e) 2 



[31] 2 [23]!- 



[31] [23] 



2-e 



[31] 2 [23] 2 - 



J(2,l -3e,e) 



2e(l-e) e(l-e) [12]e(l - e) 



J(2,2-3e,e-l) 



[31H23] 1 "' [31][23] 2 ~ e [31] 2 [23] 2 " e . 

1 - 2e N 1 - 2e 

+ 777777—77777 </(! + e, 1 - 3e, 2) - [oil2 J(e, 1 - 3e, 2) 



+ 



[12]i"«[31] 

1 - e 
[31][23] 2 " e 



[12]!- e [31] 2 ' 



+ 



l-e 



[12](l- e ) 



[31PP23] 1 "' [31] 2 [23] 2 - 



^ J(2,l-2e,0) 



e — e 



[31] 2 [23] 



— J(l,l-3e,l + e) 



(7) 



Here we use the notation 



A(a\, 012,0:3) = a(oi)a(o 2 )a(o 2 ); a(a) = 



T(D/2 - a) 
T(a) 



The second and the third integrals (5) and (6) are given in Appendix A, Eqs. (A.l) and (A. 2), 
respectively, and have a similar structure as the first integral (7). As a consequence, the total 
Integral lb can be organized into the sum 



D+2 
k=D-l 



(8) 



where 1^ represents the terms with 7(01,02,03) with a\ + o 2 + 03 = k. 

4 Integrals with the sum of indices D — 1 

The result for 1^ and for I^ D+1 ^ + I^ D+2 ^ is written in Appendix B in terms of D-dimensional 
J-integrals and then reduced to explicit expressions by using the uniqueness relation and its 
variants for the integrals. Below we consider the sum I^ D ~ l \ 



r(D-i) 



A(l,l-2e,l) 



4e(l + 2e)(l-2e) 2 \ [31] 2 [12] 



1 + 5e - 6e 2 . 

— J(e,2-3e,l) 



+ 



1 - 7e + 10e 2 1 - 4e + 5e 2 (1 - 4e + 5e 2 )[21] 



[31] 2 [23] 



l-e 



[31] [23] 



2-e 



[31] 2 [23] 2 - 



+ 



+ 



2-6e + 4e 2 2 - 4e + 2e 2 [12](2 - 4e + 2e 2 ) 



[31] 2 [23]!- 
4e - 4e 2 



[31] [23] 



2-e 



[31] 2 [23] 2 - 



J(l,2 - 3e,e) 
J(2,l -3e,e) 



[31H23] 1 -' [31][23] 2 



2e-2e 2 [12](-2e + 2e 2 ) 



[31] 2 [23] 



2-e 



J(2,2 - 3e,e- 1) 
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+ 



+ 



-e + 2e 2 
+ [i2p^[31p 
1 - e - 2e 2 



J(e,l -3e,2) 



[12]!- e [31] 2 
-2e + 4e 2 



+ 



2e - 6e 2 



+ 



e-3e 2 



+ 



[12]!- e [31] 2 
-e + 3e 2 )[12]' 



J(e-l,2-3e,2) 



J(e- l,3-3e,l) 



[31H23] 1 "' [31][23] 2 ~ e [31] 2 [23] 2 - 



J(l,3-3e,e- 1) 



(9) 



The D-dimensional integrals appearing in D cannot be found by the uniqueness method or 
its variants. We will reduce them in this Section to explicit expressions and to terms proportional 
to the integrals J(l, 1,1), J(e, 2 — 3e, 1), J(l,2 — 3e, e) and specific derivatives thereof - by 
applying the integration by parts (IBP) procedure [29, 30]. Using formula (C.l) from Appendix 
C (IBP), we represent integral (9) as 

_ 1-26,1) f -l + 4e-(9/2) g 2 

6 " 4e(l + 2e)(l-2e) 2 1 [31] 2 [12]i- l ' ' j 



+ 



1 - 6e + (15/2)e 2 1 - 7/2e + 15/4e 2 (-1 + 7/2e - 15/4e 2 )[21] 



[31H23] 1 -* 



+ 



[31] [23] 



2-e 



+ 



[31] 2 [23] 2 - e 



+ 



2-6e + 4e 2 2 - 4e + 2e 2 [12](2 - 4e + 2e 2 )' 
[31] 2 [23p^ + [31][23] 2 - e [31] 2 [23] 2 - £ 

-e + 2e 2 



J(l,2-3e,e) 
J(2,l-3e,e) 
J(e,l -3e,2) 



+ 



36 - (5/2)e 2 + (3/2)6 - (5/4)e 2 + [12](-(3/2)e + (5/4) e 2 ) 



[31] 2 [23]!- e 



[31] [23] 



2-e 



[31] 2 [23] 2 " e 

1 ~ (5/2)6 2 
[12]!^[31] 2 



[12]!- e [31] 2 

J(2,2-3e,e-l) 



+ 



A(l,l -2e,l)A(l,l -3e,l + e) 
8e 2 (l + 2e)(l - 2e) 2 



-£ - (l/2)6 2 



+ 



J(e-l,2-3e,2)j 

-(6/2) - (l/4)6 2 

[12] 2 - 26 [23] 1 ~ 2e [31] ' [12] 2 - 2e [23] 2 ~ 2e [31p 

(3/2)6 + (7/4)6 2 \ 
[12] 1 ~ 2e [2S] 2 - 26 [31] J ' 

After applying formulas (C.3) and (C.4) (IBP), this formula takes the following form: 

m_D_ ^1(1,1-26,1) / -l + 56-(19/2)6 2 

b ~ 4e(l + 26)(l-2e) 2 1 [31] 2 [12]i- J[£,Z 1J 



+ 



-l + 6e- (29/2)e^ -1 + 13/2e - 41/4e 2 (1 - 13/2e + 41/4e 2 )[21] 



[31] 2 [23]! 
+ 



[31] [23] 



2-e 



[31] 2 [23] 2 " e 



J(l,2 - 3e,e) 



3e - (5/2)6 2 (3/2)6 - (5/4)e 2 [12](-(3/2)e + (5/4)e 2 ) 



[31] 2 [23] 



l-e 



[31] [23] 



2-e 



[31] 2 [23] 2 - 



J(2,2 - 3e,e- 1) 



6 



+ 



1 - (5/2)e 2 
[l^-^l] 2 

1 [12] 



[31] 2 [23] [31] [23] 2 [31] 2 [23] 2 



J(e-l,2-3e,2) 



J(l,l,l) 



A(l,l - 2e,l)A(l,l -3e,l + e) 



-2 + 7e - 6e 2 



8e 2 (l + 2e)(l - 2e) 2 
(2-4e + 2e 2 ) 



+ 



-2 + (ll/2)e-(l/4)e 2 
■e-(l/2)e 2 + _( e /2)-(l/4)e 2 



[12]- 2e [23] 2 - 2e [31] 2 + 2e [l^ 2 - 2 ^] 1 "^]^ [12] 2 - 2e [23] 2 - 2e [31] 2e 
Finally, applying relations (C.5) derived in Appendix C, we obtain the representation 



A(l,l-2e,l) 



+ 



t(D-I) __ 

b 4e(l + 2e)(l-2e) 2 \[31] 2 [12] 

-l + 9e-17e 2 -1 + 8e - (23/2)e 2 (1 - 8e + (23/2)e 2 )[21] 



5e-12e 2 T/ „ n 

— J(e,2-3e,l) 



[31] 2 [23] 



l-e 



[31] [23] 



J(l,2-3e,e) 



+ 



36 - (5/2)e 2 -(3/2)6 + (5/4) £ 2 [l_2](-(3/2)e + (5/4) e 2 



[31] 2 [23] 



l-e 



[31] [23] 



[31] 2 [23] 2 - e 

1 - (5/2)6 2 



+ 



1 



[31] 2 [23] 2 " e 

(13) M a«J(l,2-3e,e) 

(31) M 9( 3 )j(e,2-3e,l)l 



[12]i-[31] 



+ 



[12] 



A(l,l-2e,l)A(l,l-3e,H-e) 
8e 2 (l + 2e)(l-2e) 2 

2 - 4e + 2e 2 



[31] 2 [23] [31] [23] 2 [31] 2 [23] 2 
2 + 7e - 6e 2 -3 + 6e - l/2e 2 



J(l,l,l) 



+ 



[n] 1 - 2 ^} 1 - 2 ^!}^ [n] 1 - 2 ^} 2 - 2 ^!}^ 



+ 



[12]- 2e [23] 2 " 2e [31] 2 + 2e 



+ 



-4e - e 2 



[l^ 2 " 2 ^] 1 " 2 ^!]^ 



+ 



-2e - e 2 /2 



[12] 2 - 2e [23] 2 ^ 2e [31] 2e 



(10) 



5 The sum 

Summing up Eqs. (10), (B.l) and (B.2), we obtain 

A(l,l-2e,l) 



h 



5e - 12e 2 T , 

— J(e,2-3e,l) 



+ 



-l + 9e-17e 2 -1 + 8e - (23/2)e 2 (1 - 8e + (23/2)e 2 )[21] 



[31] 2 [23] 



l-e 



[31] [23] 



[31] 2 [23] 2 - 



+ 



36 - (5/2)e 2 + -(3/2)6 + (5/4)6 2 + [l_2](-(3/2)e + (5/4) e 2 



[31] 2 [23]!- 



[31] [23] 



2-£ 



4e(l + 2e)(l-2e) 2 [ [31] 2 [12] 

ill 

J(l,2 - 3e,e) 
(13) M a«J(l,2-3e,e) 
+ ,v;--.-:-(31)^ ) J(e,2-3e,l)l 



[31] 2 [23] 2 - 

1 - (5/2)e 2 



[12]!- e [31] 2 



1 f 



-19/4 



+ 



1 

~2 
1/4 



1 1 

+ 



[12] 



[31] 2 [23] [31] [23] 2 [31] 2 [23] 2 



+ 



-11/4 



+ 



-9/4 



+ 



J(l,l,l) 
5/2 



8 \[12][23][31] 2 [12][23] 2 [31] [12] 2 [23][31] [12] 2 [23] 2 [23] 2 [31] 2 [12] 2 [31] 2 



A(l,l-2e, 1)A(1, l-3e,l + e) 
8e(l + 2e)(l - 2e) 2 
-1/2 



-7/2 



+ 



[l^ 1 - 2 ^] 1 - 2 ^]^ 
3/2 



+ 



-2 



+ 



[12]l-2e[23]2-2e[ 3 l]l+26 

3/2 



+ 



[12]- 2e [23] 2 - 2e [31] 2 + 2e [12] 2 - 2 * [23] 1 ~ 2e [31] [12] 2 - 2e [23] 2 ~ 2e [31p 
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After applying formulas (C.6) which were obtained by an application of the Gegenbauer poly- 
nomial technique (GPT, cf. Ref. [35]), we obtain for the terms proportional to J(e, 2 — 3e, 1) and 
J(l, 2 — 3e, e) in It the following explicit result: 
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The explicit results for the terms I^ J ' l \ I^' 2 ^ proportional to derivatives of 7(1,2 - 
J(e, 2 — 3e, 1) in If, are written in Appendix D in Eqs. (D.l) and (D.2), respectively. 
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6 The final result for diagram (6) 

Substitution of expressions (12), (D.l) and (D.2) into the sum (11), and its expansion in terms 
of e, is performed in Appendix E, Eq. (E.2). Performing explicitly the derivatives appearing in 
Eq. (E.2), we obtain the final result for the diagram (b): 
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The result for the Davydychev integral J(l, 1, 1) is in Ref. [37]. A new integral representation 
for it has been found in Ref. [9]. The term linear in e of the J(l, 1, 1) integral in D = 4 — 2e 
dimensions was obtained in Ref. [38], while all-order e-expansion was derived in Ref. 



7 The final result for Lcc correlator 



In this section we write the total result for planar two-loop Lcc vertex. It contains contribution 
of five planar diagrams (a) — (e). The first contribution that corresponds to diagram (a) has 
been found in Ref. [9] . We re-present the final formula of Ref. [9] in the following form 
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The total contribution of diagram (a) to the Lcc correlator is this integral multiplied by a weight 
factor coming from combinatorics of Feynman rules, trivial group algebra, and Grassmanian 
nature of the ghost fields. The contributions of all the diagrams (a), (b), (c) to the effective 
action have common factor ig A N 2 /tt 8 in front of them. We concentrate on the relative weights 
caused by degrees of 1/2 and signs. This weight factor for diagram (a) is 1/2 16 . 

The final result for diagram (b) is Eq. (13). We re-present it in the following form 
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This integral should be summed with another integral obtained from Eq. (15) by exchanging X2 
and X3. Thus, the total contribution of diagram (6) to the Lcc correlator at two loop level is 
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The weight of diagram (6) in the full result is — 1/2 13 . 

The result for diagram (c) is obtained in a different way than for diagrams (a) and (b) and 
the details of calculation of the diagram (c) can be found in Ref. [59]. It takes the form 
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This integral should be summed with another integral obtained from Eq. (17) by exchanging x<i 
and X3. Thus, the total contribution of diagram (c) is 
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The weight of diagram (c) in the full result is 1/2 15 . 

Thus, the contribution of first three diagrams (a), (b) and (c) to the two- loop effective action 
can be presented in the following form 
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Diagrams (d) and (e) are easy to calculate since they contain propagator-like insertions only. 
They can produce logarithms only (due to the shift of the indices in the propagators) and 
cannot produce Davydychev integral J(l, 1, 1). These two diagrams are divergent separately in 
the ultraviolet region, however in the maximally supersymmetric Yang-Mills theory their poles 
cancel each other. Thus, the contribution of diagrams (d) and (e) to the two-loop effective action 
can be presented in the following form 
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Thus, the full result for the two-loop Lcc correlator in the maximally super symmetric Yang- 
Mills theory in four spacetime dimensions is a combination of V^ a \ V {b \ V® and V {d+e ^ with 
the correpsonding weights and has the form 



/ • 4 at2 
d i x 1 d i x 2 d 4 x 3 j^f abc L a (x 1 )c b (x 2 )c c (x 3 ) V^( Xll x 2 ,x 3 ) 



where 



29/3 



,22,23) = V ( - a+b+c \xi, 22,23) + y( d+e ) (21, 2 2 ,2 3 ) = 

32/3 29/3 -52/3 _ -40/3 _ -52/3 



+ 



+ 
+ 



+ 



+ 



[12] 2 [23] 2 [12] 2 [31] 2 [23] 2 [31] 2 [12][23][31] 2 [12][23] 2 [31] [12] 2 [23][31] 

1 1-13 -13 -2[12] 



+ 



+ 



+ 



+ 



[12] 2 [23] [23] [31] 2 [12] [31] 2 [12] 2 [31] [23] 2 [31] 2 



-2[31] 



+ 



+ 



+ 



12 [23] 
[12] 2 [23] 2 ' [12] 2 [31] 2 
6 18 



+ 



+ 



J(l,l,l) 
10 



[12] 2 [23] 2 [12] 2 [31] 2 [23] 2 [31] 2 [12] [23] [31] 2 
-14 6 



+ 



+ 



+ 



-16 



+ 



[12][23] 2 [31] [12] 2 [23][31] 
12 -16 



ln[12] 



+ 



+ 



16 



[12] 2 [23] 2 [12] 2 [31] 2 [23] 2 [31] 2 [12] [23] [31] 2 
28 -16 

+ ^77^ + 



18 



+ 



[12][23] 2 [31] [12] 2 [23][31] 
-6 -2 



In [23] 



+ 



+ 



[12] 2 [23] 2 [12] 2 [31] 2 [23] 2 [31] 2 [12] [23] [31] 2 
.„'!.,,. + ttt^LtttI ln[31] 



[12][23] 2 [31] [12] 2 [23][31] 



13 



For the reference, the one-loop contribution to the effective action is 

J d 4 x 1 d 4 x 2 d 4 x 3 l -^f abc L a (x 1 )c b (x 2 )c c (x 3 ) vW( Xl ,x 2 ,x 3 ), 

where 
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8 Conclusions 

In this paper we presented the calculation of the second diagram [diagram (&)] of the five planar 
two-loop diagrams to the Lcc vertex depicted in Fig. 1 for a general Yang-Mills theory. The 
Lorentz structure of the corresponding integrand is complicated in comparison with the integrand 
of the previously calculated diagram (a) [9]. In order to apply the same methods as in the 
diagram (a), we first performed in the integrands the Lorentz algebra for certain subproducts 
in the numerators, then performed one of the two D-dimensional integrations, then finished 
the Lorentz algebra, and finally performed the second D-dimensional integration. It is a long 
procedure that requires certain computer resources. In principle, it is possible to use another 
trick to reduce the calculation [59], a trick which we applied in the calculation of diagram (c) 
[59]. Neither the result for the diagram (b) obtained here, nor the result for the diagram (a) 
obtained in our previous work [9] depend on any scale, infrared or ultraviolet, in complete 
correspondence with the naive index counting arguments, taking into account that none of three 
diagrams (a, b, c) contains any divergent subgraphs due to transversality of the gauge propagator 
in the Landau gauge. From the point of view of R-operation theory the scale independence of 
this vertex in the Landau gauge in all orders of perturbation theory has been explained in Refs. 
[1, 2, 3]. By using the ST identity, other (leading- iV) two-loop vertices can be derived from 
the (leading-TV) two-loop Lcc vertex once the latter is known. For example, it will be possible 
to derive the four-point off-shell correlator, and consequently reproduce the known result for 
the four-point gluon amplitude [18, 19] and the anomalous dimensions of twist-two operators 
[40, 41, 42, 43, 44]. 

It is possible to look at these results from different points of view. On the one hand, singular 
parts of the diagrams (d) and (e) are proportional to the one-loop result for the Lcc correlator 
written in Ref. [9]. Moreover, the transversality of the gauge propagator must be conserved by 
the radiative correction, since it is a well-known fact that the gauge fixing term does not obtain 
any quantum correction [15]. This means that the sum of singular parts of diagrams (d) and 
(e) is proportional to the one-loop contribution, with a coefficient that is singular in a general 
nonsupersymmetric massless gauge theory, but in this particular theory it is a finite factor since 
the poles in e cancel in the sum of the diagrams (d) and (e) due to the Af = 4 supersymmetry. 
This factor is logarithm of a ratio of the spacetime intervals. This factor appears due to shift 
of indices by multiplies of e in the propagator of the fields. This shift is typical for massless 
theories in MS scheme. In the nonsupersymmetric case, that singular number can be absorbed 
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into the gauge coupling to organize the bare coupling. Then, the bare coupling together with 
the logarithm of ratio of the distance to the scale, leads to the running (effective) coupling. It 
means that in D dimensions the massless nonsuper symmetric gauge theory is a conformal gauge 
theory in terms of the running effective coupling (formed from the bare coupling) and dressed 
mean fields. In other words, we think that, in terms of the running couplings and dressed mean 
fields, the arguments of Ref. [1] can be applied without modifications to the massless QCD. The 
role of the renormalization group (RG) scale could be the coordinate of the moduli space of the 
theory. 

On the other hand, QCD is asymptotically free theory and for short distances we can consider 
it as a theory with zero beta-function (running of the effective charge is almost absent) and the 
conformal structure could be restored at short distances by the method of conformal theory 
in terms of dressed mean fields. Further, it would be interesting to investigate the relation of 
off-shell correlators of dressed mean fields and correlators of instantons in multi-dimensional 
theories. To introduce masses in the theory, we need to use softly broken super symmetry, in 
which the couplings are spacetime-independent background superfields. The relation between 
the RG functions of softly broken and rigid theories was found out in Refs. [45, 46, 47]. The 
relation between the correlators of softly broken and rigid theories can be found by a trick of 
general change of variables in superspace [48] . 

To restore the conformal structure of all the effective action, we need to solve the ST identity. 
Algorithm for solving the ST identity could be applied to various theories, such as the N = 8 
supergravity [49, 50, 51, 52], Chern-Simons theory near the RG fixed points [53], massless gauge 
theory near fixed points in the coupling space, topological field theories in higher dimensions, 
finite M = 1 supersymmetric theories [54, 55, 56, 57, 58]. We further note that, in addition 
to the five planar two- loop diagrams of Fig. 1, there is one nonplanar diagram (the nonplanar 
variant of the diagram (a)) which is suppressed in the planar limit of the large 't Hooft coupling 
by the factor l/N. 
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Note added 

After publication of the first version of this paper the "dual conformal symmetry" of the mo- 
mentum space (conformal symmetry of four-point correlator of gluon in the momentum space) 
found in Ref. [20] has been exploited further in Ref. [60] up to four loop level to detect all 
off-shell correlators that will contribute to four-point amplitude. That symmetry appears also 
on the strong coupling side as it has been found by Alday and Maldacena in Ref. [61] via 
AdS/CFT correspondence. As we have already mentioned in Introduction, Fourier transform 
of the UD integrals of the momentum space can be expressed in terms of UD integrals in the 
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position space. This observation could be a bridge between the dual conformal symmetry of the 
momentum space and conformal symmetry of the effective action of dressed mean fields in the 
position space. 



Appendix A 



In this Appendix we present the result for the second and the third integrals (5) and 
of D-dimensional integrals J. The result for the second integral is 



h,2 = 



A(l,l-2e,l) 



20e)(l-e 2 



8e(l-2e) 2 ( [31][12]- e 
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The third double integral (6) is simpler, the integration over y can be taken immediately, 
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i + lzif! | (^-2 £ 2 )[23] l-2 e -2 e 2 \ J(1 + e2 _ 3£l 
[12]^[31] 2 [12] 1 ~ e [31] 2 [12]!- e [31] J J ^ e ' z ^ 
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+ ftttti , r „, 10 2 - 3e, 1) - Mn11 _ froi19 J(l + e, 2 - 3e, 



[l^-^l] 2 

2e - 2e 2 

J(l + e, 3 - 3e, 0) + J(e, 3 - 3e, 



+ 



[12]!^[31] 2 
26 - 2e 2 
" [12]-[31] 2 

' (2e-2e 2 )[23] 

[12]" e [31] 2 [12]" e [31] [12]- £ - 1 [31] 2 / 
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(A.2) 



In Eqs. (A.l) and (A.2), we do not write single integrals with one of the indices zero and sum 
of other two indices equal to D, D + 1 or D + 2. Such terms are proportional to delta-functions 
in position space multiplied by pole in e and will disappear in the final expression (13). 



Appendix B 

In this Appendix we write the results for the terms contributing to the diagram (b) with J's 
whose sum of indices is k = D, D + 1, D + 2. These terms appear in 1^ which is the sum of Eqs. 
(7), (A.l) and (A.2). The integrals in are calculated by using the uniqueness method 

(D) _ A(M-2M) f / 10 - 186 - 126 2 3-4e-10e 2 (2e - 2e 2 )[23] \ 

6 4e(l + 2e)(l-2e) 2 1 I [31] 2 [12]^ ^ [31H12] 1 "' [31] 2 [12]!^ J^ 1 ^ 6 '^ oe >^ 



18 



+ 



-3 + 3e + 12e z 



[31] [12] 



l-e 



J(2 + e,l-3e,l) + 



2e - 4e 



[31] 2 [12] 



— J(2 + e, 1 -3e,0) 



-4e + 8e 2 , . -2e + 2e 2 , 

f t^ttt^ ^(1 + e, 2 - 3e, 0) + ^^-,7(2, 1 - 3e, 1 + e) 



[31] 2 [12] ] 



+ 



+ 16e - 6e 2 -e + 3e 2 (-2e + 5e 2 ) [21] 
- + : + 



[31] 2 [23]" e [31][23]!- e 
+ 



[31] 2 [23] 



l-e 



J(l,2-3e,l + e) 



(7e - 14e 2 )[21] (e-3e 2 )[21] [12] 2 (-e + 3e 2 



[si] 2 ^] 1 -^ 



[31] [23] 



2-e 



8 - 20e + 12e 2 -2e + 4e 2 



[31] 2 [23]- 



[31] [23] 



l-e 



[31] 2 [23] 2 " e 
J(2,2-3e,e) 



+ MJ^r^ J(l + e, 1 - 3e, 2) + 



[12]i"«[31] 
2e - 4e 2 



/ -2 + 2e t -2 + 2e [12](2 - 2e) 
I __ . h r „, lor „„ n _ e + 



V[31][23] 2 - e [31] 2 [23] 



[31] 2 [23] 



2-e 



+ 



[12]!^[31] 2 



J(e,3 - 3e,0) + 



-2 + 6e 2 



+ 



(-6) [23] 



+ 



-2 + 5e - 2e 2 



+ 



+ 



[l^-^l] [12]!-[31] 2 [12]" e [31] 2 
(-2e + 2e 2 )[23] e - 2e 2 4-lle + 2e 2> 



+ 



+ 



v [12]!- e [31] 2 [12] 1 ~ e [31] [12]~ e [31] 2 
'4-18e + 26e 2 -e + 4e 2 (3e - 10e 2 )[21] ' 



[31] 2 [23]- 



+ 



[31][23]i- 



+ 



+ 



[31] 2 [23]!- e 
[12]!- e [31] 2 



J(2,l-2e,0) 
J(e,2 - 3e,2) 

J(e, 3 — 3e, 1) 

J(l,3 - 3e,e) 



J(e-l,3-3e,2) 



+ 



(-2e + 6e 2 )[21] (-e + 3e 2 )[21] [12] 2 (e-3e 2 



[31] 2 [23] 



l-e 



+ 



[31] [23] 



2-e 



[31] 2 [23] 2 - 



J(2,3-3e,e-l) 



e - 3e 2 -e + 3e 2 

+ 77^7^7737 «/(0, 2 - 36, 1 + e) + — — — - J(0, 3 - 3e, e) 



+ 



[31H23] 1 "' v ' ' ; [31] 2 [23] 

-l + 3e-2e 2 -1 + 3e - 2e 2 (1 - 3e + 2e 2 )[12] ' 



+ 



+ 



[31] [23] 2 ^ [31H23] 1 " 6 [31] 2 [23] 2 ~ 

A(l,l-2e,l)A(l + e,l-3e,l) f 



J(l,2 - 2e,0) 
-2e + 10e 2 



8e 2 (l + 2e)(l - 2e) 2 



3 + 7e + 3e 2 



+ 



2 - 4e - 5e 2 



[12] i- 2 * pa] 2 - 26 [3i] x + 2e 



+ 



[ 12 ]2-2e[23]l-2e[ 13 ]l+2e 

2-4e 

2e [OQl2-2e[Q112+2e 



[12]- 2e [23] 2 - 2e [31] 



+ 



+ 



-5 + 8e - 18e 2 



[12] 2 ~ 2e [23]- 2e [13] 2 + 2e 



+ 



-6e + e 2 



e - (l/2)e 2 



+ 



[ 12 ]2-2e[23]2-2e[ 31 ]2e [31]2[23] 



[31]2[ 12 ]2-2e 
2rool2-2e [ 



A 2 (l,l-2e,l) 



1 - e 



+ 



l-e 



l-e 



4e 2 (l + 2e)(l-2e) 2 I [12] 1 - fi [23] 2 - e [31] [l^ 1 -^] 1 -^] 2 [12]- e [23] 2 ~ e [31] 2 

The terms in the sums and /^■ D+2 - ) can be calculated also by using the uniqueness method, 

but only after having represented the J-terms as derivatives of integrals J(oi\, 02,03) with 
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In this Appendix we provide formulas used in Sec. 4. The following formula results from applying 
the corresponding IBP procedure: 
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+ (e-l) J(l,2-3e,e) 
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[12] 2 - 2e [23]- e [31]- 1 + 2e , 

Other useful formulas can be taken from Ref. [9] (Appendix A there, obtained by IBP procedure): 
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The following useful relation can be directly derived: 
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2e(l-3e) [23] 2 " 2e [12]- e [31]- 1+2£ 
For J(2 — 3e, 1, e) we use specific formulas of Ref. [9] which were obtained by applying the GPT: 
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Appendix D 

In this appendix we apply formulas (C.6)and obtain for the term in Eq. (11) proportional to 
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Appendix E 



Collecting (12), (D.l) and (D.2) in Eq. (11), we obtain the final expression for the diagram (b) 
in D dimensions 
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Taking the limit e — ► 0, the poles in e disappear and the result is 
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